Description of a quantum convolutional code 
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We describe a quantum error correction scheme aimed at protecting a flow of quantum information 
over long distance communication. It is largely inspired by the theory of classical convolutional 
codes which are used in similar circumstances in classical communication. The particular example 
shown here uses the stabilizer formalism, which provides an explicit encoding circuit. An associated 
error estimation algorithm is given explicitly and shown to provide the most likely error over any 
memoryless quantum channel, while its complexity grows only linearly with the number of encoded 
qubits. 

PACS numbers: 03.67.Pp, 03.67.Hk, 03.67.-a 



In recent years, the discovery and development of 
quantum computation and communication has shed new 
light on quantum physics. The potential applications 
of these new fields encompass a wide variety of sub- 
jects, ranging from unconditionally secure secret key gen- 
eration protocols to efficient integer factoring algo- 
rithms [2} or enhancement of communication complex- 
ity 0. However, the practical realization of such proto- 
cols and algorithms remains a very involved task mainly 
because of the inherent instability of quantum superpo- 
sitions |4| as well as intrinsic imprecisions of the physical 
devices that process quantum information. These errors 
wipe out the quantum superpositions together with en- 
tanglement, which are usually seen as key resources of the 
power of quantum algorithms and protocols Hence, 
protecting the quantum nature of information became 
one of the most important challenges to prove the feasi- 
bility of quantum computers. The discovery of quantum 
error correction schemes 0, notably opened the future 
of large scale quantum information processing: a certain, 
but unfortunately very small, degree of imprecision can 
be tolerated at each step of a quantum transformation 
and still allow a speed-up over classical information pro- 
cessing However, building a fault-tolerant quan- 
tum computer remains largely out of reach of the present 
day practical realizations, principally because of the large 
number of physical qubits required to account for the er- 
ror correction. 

On the other hand, quantum cryptography and more 
generally the field of quantum communication seems 
more promising in a near future. Some quantum key dis- 
tribution protocols have been implemented and the asso- 
ciated devices seem to be close to commercialization ^(j . 
Within this context, we construct a new family of codes 
- quantum convolutional codes — aimed at protecting 
a stream of quantum information in a long distance com- 
munication. They are the correct generalization to the 
quantum domain of their classical analogues, and hence 
inherit their most important properties. First, they have 
a maximum likelihood error estimation algorithm for all 
memoryless channels with a complexity growing linearly 
with the number of encoded qubits. This is an impor- 



tant issue since finding the most likely error — a strat- 
egy which allows to determine the most likely sent code- 
word — is in general a hard task: for a generic family of 
block codes with constant rate, the maximum likelihood 
error estimation algorithm has a complexity growing ex- 
ponentially with the number of encoded qubits. Hence, 
generic block codes rapidly require to employ suboptimal 
error estimation procedures which, as a consequence, do 
not exploit the whole error correcting capabilities of the 
code. Moreover, our algorithm can easily handle vari- 
ations in the properties of the communication channel 
(i.e. a change in the single qubit error probabilities). 
The second advantage of quantum convolutional codes is 
their ability to perform the encoding of the qubits online 
(i.e. as they arrive in the encoder). Thus, it is not neces- 
sary to wait for all the qubits to be ready to start sending 
the encoded state through the communication channel: it 
reduces the overall processing time of the qubits which 
is an additional source of decoherence. Note that an at- 
tempt at defining quantum convolutional codes has been 
made some time ago [ill 1 1 2| . but missed some crucial 
points concerning the error estimation algorithm as well 
as error propagation properties. 

In this letter, we deal with a specific example drawn 
from our general theory. We construct a quantum con- 
volutional code achieving a rate equal to 1/5: we explain 
how to encode and decode a stream of qubits efficiently, 
and we expose the maximum likelihood error estimation 
algorithm. This will give all the necessary intuition to un- 
derstand how to generalize the present results to a wider 
framework 13]. 

Description of the code - The particular code we 
wish to present is best described by using the stabilizer 
formalism |14j . This provides a simple way to under- 
stand the encoding and decoding operations. Moreover, 
the error syndromes can be easily identified, which con- 
siderably simplifies the description of the error estimation 
algorithm. We use the following standard notations for 
the Pauli operators acting on a single qubit: 
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so that ZX = Y . The identity matrix will be denoted by 
/. Since convolutional codes are designed to deal with a 
stream of information qubits, the number of generators 
of the stabilizer group will possibly be infinite. However 
in practice, transmission starts and ends at a given time, 
which means that we only consider generators made of a 
finite number of Pauli operators. 

The code subspace is described by the generators of its 
stabilizer group, S. These generators are given by: 



M 


= XZ I I I I I I 


Mi 


= Z XXZ I III..., 


M 2 


= I Z XXZ I II..., 


M 3 


= 11 Z XXZ I I 


M 4 


= 111 Z XXZ I 


M 4i+j 


= I® 5t ®M , 0<i, 1 <j < 4 , 


M x 


= ... I I I I Z X. 



It is easy to check that all the generators commute and 
are independent. Thus, the code subspace (i.e. the largest 
common eigenspace of the generators with eigenvalue +1) 
is non-trivial. 

An important point to address when considering sta- 
bilizer codes is the ability to manipulate encoded infor- 
mation. Namely, we want to find the encoded Pauli op- 
erators Xi, Zi corresponding to logical qubit i. These 
operators must satisfy the following relations: 



X u Z, G N(S) - S, (3) 

Ji:Xj. = 0, (4) 

[Zi,Zj} = 0, (5) 

[X l ,Z 3 } = 0, t^j, (6) 



where N(S) denotes the normalizer of S. Equation 
states that the encoded Pauli operators leave the code 
subspace globally invariant, but have a non-trivial action 
on its elements, while the Equations (1 lKill ensure that 
manipulating qubit i does not affect other qubits. There 
exists a great choice of different sets of such operators, 
however they are not equivalent in the perspective of ef- 
fectively manipulating the encoded quantum information 
in an easy way: in practice only those with a small num- 
ber of terms different from the identity are useful. For 
our particular example, such set exists and has a struc- 
ture invariant by a shift of five qubits: 

X X =IZIXIZII..., 
~Zx = I Z Z Z Z Z I I , . 
~X n = J® 5 "® X l7 n> 1, ( ' 

Z n = /® 5n (g)Zi, n>l. 

Hence, a unitary transformation on a single encoded 
qubit will in general be implemented by a unitary trans- 
formation on five physical qubits. 

At this point, one can wonder what in this code dif- 
fers from a generic block code. The answer to this ques- 
tion comes from the particular structure of the stabilizer 



generators: beside M and M^, the generators of the 
stabilizer group can be casted into sets of constant size 
(e.g. four), each set acting on a fixed number (e.g. seven) 
of consecutive qubits. In addition, each set has a fixed 
overlap (e.g. of two qubits) with the set immediately be- 
fore and immediately after. This very peculiar struc- 
ture defines quantum convolutional codes and we can 
prove [lii ] that this implies the possibility of online en- 
coding and the existence of an efficient error estimation 
algorithm. 

Encoding circuit - As explained in D. Gottesman's 
Ph.D. thesis ^ij, there are various ways to realize the 
encoding into the code subspace. However, for convolu- 
tional codes, they are not equivalent: standard encod- 
ing circuits usually require to wait until the last 'to-be- 
protected' qubit has been obtained before sending the 
encoded state. In this section, we explain how to take 
advantage of the structure of the stabilizer generators to 
overcome this limitation and encode the qubits online. 
We first exhibit a map from the computational basis of 
the 'to-be-protected' qubits to a basis of the code sub- 
space. As a second step, we derive the quantum circuit 
implementing this map in a unitary way. 

More precisely, consider the following set of states: 

{|?A(ci,C2,C3,...))} Ci6{01} = (8) 

{p\o, o, o, o, o, a, o, o, o, o, C2 , o, o, o, o, c 3 , . . .)} Cie{0 , 1} , 

where P = Yl^I + Mj)/y/2 is the projection operator 
onto the code subspace. Since Zi commutes with all the 
generators of the stabilizer group, the following equation 
holds for any element of the set: 

Z l P\0, 0, 0, 0, 0, a, 0, 0, 0, 0, c 2 , 0, 0, 0, 0, c 3 , . . .) = (9) 
(-1) C 'P|0, 0, 0, 0, 0, ci, 0, 0, 0, 0, c 2 , 0, 0, 0, 0, c 3 , . . .) . 

This implies that {\ip(c\, c 2 , c 3 , . . -))}c 4 e{o,i} is an or ~ 
thonormal basis of the code subspace. Hence, the natural 
encoding consists in mapping the computational basis of 
the 'to-be-protected' qubits, {|ci,c 2 ,c 3 , . . .)} Ci e{o,i} 
the basis {\ip(c 1} c 2 , c 3 , . . .))} Ci e{os}- 

In practice, to encode a stream of qubits qi, we first 
add to it ancillary qubits in the |0) state such that the 
'to-be-protected' qubit i is now at the position 5i + 1. 
Then, we need to implement P for these specific input 
states as a unitary transformation onto the whole Hilbert 
space. This can be done in full generality as explained 
in 14], and gives the encoding circuit of Fig. ^ From 
this simple example, it is easy to understand that the 
possibility of online encoding for quantum convolutional 
codes is a consequence of the finite extension of the sup- 
port of the generators of the stabilizer group and of the 
encoded Pauli operators. Also note that alternative en- 
coding methods can be found and can be relevant when 
considering some specific applications, but these issues 
are beyond the scope of this letter. 
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Error propagation and online decoding - Due to 
their very specific nature, convolutional codes propagate 
information contained in a given qubit to its successors 
(see again FiG.^]). During the decoding process (i.e. the 
inverse of encoding) this can actually become a problem: 
an error affecting a finite number of qubits before de- 
coding can propagate through the decoding circuit and 
finally affect an infinite number of qubits. Such error 
is called catastrophic. It is worth mentioning that this 
issue is to the quantum domain: classical convolutional 
encoders might have catastrophic errors 0, For- 
tunately, in both cases, non-catastrophic encoders exist. 
More precisely, given a specific encoder one can find a 
procedure to determine whether it has catastrophic er- 
rors or not. For classical codes this is a well known re- 
sult established by Massey and Sain 01 ■ For the quan- 
tum domain, we can show 13] that the circuit of non- 
catastrophic encoders fulfills the following requirement: 
its gates form a finite number of layers and commute with 
each other inside a layer. The idea behind this theorem 
is simple. In general, an error affecting some qubits will 
propagate to all the other qubits involved in a gate with 
the erroneous ones. When those qubits are further used 
in other gates the error continues to propagate until no 
more gates are applied. The commutation relation to- 
gether with the finitencss of the number of layers ensures 
that, for any finite size error, only a finite number of gates 
will enter in the propagation process. Thus all errors are 
non-catastrophic. Fig. [5] illustrates this 'pearl-necklace' 
structure for our example, and thus proves that our rate 
1/5 quantum convolutional code is non-catastrophic. 

Moreover, it can be shown that this condition implies 
the existence of a forward decoding scheme: there is no 
need to wait for the last qubit to start decoding (see 
FiG.EJ). For non-catastrophic codes, both encoding and 
decoding can be done online |l3|. 

Maximum likelihood error estimation - An error 
correcting code aims at protecting information sent over 
a noisy communication channel by letting the receiver 
infer which error possibly affected the information. This 
is the role of the error estimation algorithm. On aver- 
age, the correct information is most often retrieved when 
the estimated error coincides with the most likely error. 
Thus, it is both of theoretical and practical relevance to 
have an efficient maximum likelihood error estimation al- 
gorithm for our quantum convolutional codes. In this sec- 
tion, we exhibit such algorithm. It is indeed the quantum 
analogue of the well-known Viterbi algorithm for classi- 
cal convolutional codes. The Viterbi algorithm realizes 
a maximum likelihood error estimation on all memory- 
less channels with a complexity linear in the number of 
encoded bits. This explains why classical convolutional 
codes are so widely used for reducing the noise on com- 
munication channels. 

Our algorithm for quantum convolutional codes pro- 
cesses the information obtained through the syndrome in 



order to infer the most likely error. The circuit for ob- 
taining the syndromes follows the usual phase estimation 
scheme: an ancillary qubit is prepared in the |0) state; 
undergoes a Hadamard transform; controls the applica- 
tion of one of the generator Mj of the stabilizer group; 
again undergoes a Hadamard transform; and is finally 
measured in the {|0), |1)} basis. Then, the algorithm up- 
dates a list of maximum likelihood error candidates by 
looking at a small number of syndromes at a time, and 
by taking local decisions. It is preceded and followed by 
appropriate initialization and termination steps. 

The initialization step lists all error candidates, {Ej}j, 
for the first two qubits which are compatible with the syn- 
drome Mq. There are exactly 8 = 4 2 /2 of them (there 
are 4 2 different operators with support on the first two 
qubits, but the constraint associated with Mq divides this 
set into two equal parts). This list constitutes the input 
of the main loop of the algorithm. At step i, the algo- 
rithm constructs a list of some most likely error candi- 
dates, {Efij compatible with the syndromes Mq to Mn. 
Each candidate Ej is thus specified only on qubits 1 to 
hi + 2. The crucial point of the algorithm is to maintain a 
fixed size of this list, and hence to avoid the exponential 
blow up that would arise when listing all error candi- 
dates compatible with these syndromes. More precisely, 
Ej is a most likely candidate whose restriction on qubit 
5i+ 1 and hi + 2 is prescribed by the index j running over 
the set of 16 possible errors affecting those two qubits. 
The computation of any error candidate E^ is easily 
achieved provided {Ej}j: consider the set of all possible 
extensions of the error candidates E* to qubit hi + 3 to 
hi + 7 with the prescribed error k at position hi + 6 and 
hi + 7. It is easy to check that any such element is now 
compatible with syndromes M to M^ i+1 y The specific 
candidate E]^ 1 is chosen to be the most likely operator 
among the elements of the latter set (in case of tie, one 
is chosen at random). This procedure is continued until 
reaching M^, which again selects half of the candidates. 
The termination of the algorithm outputs the most likely 
candidate among the remaining ones. This constitutes 
the most likely error given the value of all the syndromes 
for the received stream of qubits 

The main property used to prove this fact is related 
to the structure of the generators of the stabilizer group: 
the value of the syndromes associated to Ma+i to 
depend on the syndromes Mq to only through the 
error operators at position hi + 1 and hi + 2. Thus, taking 
a sequence of local decisions allows to construct a list of 
error candidates among which one will coincide with the 
most likely error until qubit hi + 2 while maintaining a 
linear complexity of the algorithm as the number of en- 
coded qubits increases. Note that, the error maximizing 
the likelihood is known when the last syndrome is mea- 
sured. Hence, it is in principle necessary to wait till the 
end of the transmission to actually correct the estimated 
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error. However, as for the classical Viterbi algorithm, nu- 
merical simulations show that the different candidates at 
a given step coincide with the most likely error except on 
their last few positions. Thus, in practice it is possible to 
estimate the error online. In addition, we want to stress, 
that without increasing its complexity, this algorithm can 
take into account all memoryless quantum channels even 
if the single qubit error probabilities are not constant in 
time. For example, one could imagine that the qubits 
are photons sent through an optical fiber, and that the 
probabilities are evaluated by sending probe photons con- 
taining no useful information. Finally, as the codes de- 
scribed here are the exact translation to the quantum 
setting of the classical convolutional codes, one can also 
derive suboptimal error estimation algorithms (for their 
classical analogues see |0, 0] ) . Most importantly, quan- 
tum convolutional codes can be decoded iteratively and 
allow quantum turbo decoding [l3| . 

Conclusion - In this letter, we presented the the- 
ory of quantum convolutional codes by an example. We 
gave explicitly the associated encoding and decoding cir- 
cuits, as well as a low complexity maximum likelihood 
error estimation algorithm. We believe that such codes 
could be used to reduce errors for long distance quantum 
communications provided that we are able to perform a 
small and fixed number quantum gates with good fidelity. 
Moreover, the tools developed for quantum convolutional 
codes can be used to translate other families of classical 
codes to the quantum domain, like for instance low den- 
sity parity check codes. 

Part of this work was done when H.O. was visiting 
the Perimeter Institute and the Institute for Quantum 
Computing in Waterloo, Canada. Useful discussions with 
J. Kempe, R. Laflamme and D. Poulin are gratefully ac- 
knowledged. 
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FIG. 1: Beginning of the circuit realizing the encoding once 
the ancillary qubits have been added to the stream containing 
the initial quantum information (qubits qi,qz, ■ ■ .). H is the 
Hadamard transform, and the dot represents the control qubit 
for a given gate. The circuit is run from left to right. When 
all the transformations have been performed for a given qubit, 
it can be sent through the communication channel. 
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FIG. 2: Left: The encoding circuit of Fig.Qwhere consecutive 
blocks of operations have been placed in different orders and 
the appropriate commutators introduced. There are 6 layers 
of gates in this circuit and in each layer all the gates com- 
mute with each other. It is was we call the 'pearl-necklace' 
structure. Right: Corresponding decoding circuit obtained by 
running the modified encoding circuit backward. In this form 
it is obvious that the decoding circuit has a structure allowing 
a forward decoding. 



